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1 Introduction 

A well known consequence of the fundamental equation of rotational dynamics is that 
the period of small oscillations of a physical pendulum is given by 




T= — = 27rW^ (1) 



where m is total mass of the pendulum, I its moment of inertia with respect to the 
center of oscillation O and h the distance of the center of mass from O. Then a physical 
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pendulum oscillates like a simple pendulum of length 

which is called the equivalent length of our physical pendulum. 

By the Huygens-Steiner theorem (also known as the "parallel axis theorem") it is possible 
to write 

/ = mh"^ + lo 

where Iq is the moment of inertia with respect to the center of mass. By squaring equation 
(1) we get the following quadratic relation 

h^-lh + - = (3) 
m 

When P — 4:Io/m > that equation admits two real solutions hi, /i2 such that 

hi + h2 = l (4) 

In 1817 Captain H.Kater thought to use this last relation to empirically check the 
Huygens-Steiner theorem. At this purpose he constructed his reversible pendulum con- 
sisting of a plated steel bar equipped with two weights, one of which can be moved along 
the bar. This pendulum is reversible because it can oscillate about two different suspen- 
sion points realized by two knife edges symmetrically located on the bar. By adjusting 
the movable weight, it is possible to obtain a pendulum mass configuration such that the 
periods about the two pivots coincide, the equivalent length I is the distance between the 
two knife edges and condition (4) is satisfied. 

The measurement of such a common period T, of the total mass m and of the distance I 
between the two knife edges, gives then an easy way to perform an empirical measurement 
of the earth's (apparent) gravitational acceleration g by applying formula (2). This is 
why the Kater reversible pendulum is one of the favourite instrument to measure g in 
student labs. 

Anyway there is a subtle point in this procedure which is the determination of the right 
mass configuration of the pendulum. This problem gives rise to the following two ques- 
tions: 

(1) how many possible positions of the movable weight determine a "good" mass con- 
figuration for which the periods of small oscillations about the two pivots coincide? 

(2) when a good mass configuration is realized, is the equivalent length I necessarily 
represented by the distance between the pivots? 

If the answer to the second question is assumed to be "yes" then the quadratic equation 
(3) gives precisely two possible good mass configurations since h depends linearly on the 
position X of the movable weight. These mass configurations can then be empirically 
obtained by the following standard procedure [1]: 

• by varying the movable mass position x collect two series of data (x, T), one for each 
pivot. 
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• make a parabolic fitting of the data by means of two parabolas of the following type: 



• these parabolas meet in at most two points {xi,T), {x2,T): positions Xi and X2 
determine the two desired good mass configurations. 
Such a parabolic fitting is justified by two considerations. The first one is that we are 
looking for two good mass configurations, then the fitting curves have to admit at most 
two intersection points. The second one is the empiric observation of the data which 
apparently seem to be arranged just along two convex parabolas with vertical axis. 

This is what is usually done although the right answer to the second question should 
be "no", as was firstly pointed by Shcdd and Birchby in 1907 [2, 3, 4]. Their remark seems 
to have escaped general attention, perhaps due to the fact that, if the pendulum is well 
assembled, the previous parabolas meet at points whose abscissas give almost exactly the 
good mass configurations having the distance between pivots as equivalent length. The 
latter is much easier determined than any other equivalent length associated with further 
good mass configurations of the pendulum [5] ! But what does mean "well assembled" ? 

To fix ideas consider an "ideal" Katcr pendulum consisting of an idealised masslcss 
rigid rod (x-axis) supporting two identical point masses, rrif fixed at —a and rrim at a 
variable position x. The assembly has two distinct suspension points for the oscillations 
positioned at —d/2 and at +d/2 , a sketch of this "ideal" pendulum is reported in Figure 1. 
Then d is the distance between the pivots, the center of mass is located at 



T — ax^ + bx + c 



(5) 



movable mass I 
fixed knife edgef 



fixed knife ddge 




9 



Fig. 1 Front view of idealized pendulum with two point masses and a massless rigid rod. 
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and the moment of inertia about the center of mass at b is given by 



Iq = {b + mf+{b- x) 



where m — 2m f — 2mm is the total mass of the pendulum. The moment of inertia /i 
and I2 with respect to the two pivots are 



and 



h = {^rm+i^-brm . 



(6) 
(7) 



When X determines a good mass configuration the resulting periods Ti and T2 of small 
oscillations about the two pivots, respectively, have equal values. Equation (1) gives 



\ 



m(6+f)2 


+ f{x + af 


mg 






m(6- 1)2 


+ f (^ + a)' 


mg 


b-^ 
^ 2 





= 27r, 



27r, 



{x — a + df + {x + of 


2^ 


X — a + d\ 


{x — a - 


- df + (x + af 


2^ 


X — a — d\ 



Then — 7f gives a cubic equation on the variable x. If it is assumed that 

ix -af -d'^ <Q 



(8) 



which occurs, for instance, when suspension points are the end points of the pendulum 
bar, one finds that 



{x - a) [2 (x^ + a^) - d'^ 







Its solutions are then given by 



X — a, 



X — ±\ 



-a2 



(9) 

(10) 
(11) 



which represent all the possible positions of the movable weight giving a good mass con- 
figuration for the ideal Kater pendulum. The first solution, x = a, always exists. Further- 
more, if d/y/2 > a there are two additional positions which are symmetric with respect 
to the origin i.e. the middle point of the massless bar. Recall formula (2) to obtain the 
associated equivalent lengths. For the last two symmetric solutions it gives 

/ = d 



But the equivalent length associated with the first solution is 

d 

'=2+^7 



(12) 
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which in general does not coincide with the distance d between the two pivots. 
On the other hand if (8) is not assumed and we are in the more "pathological" case of a 
pendulum such that {x — a)^ — (P > then = 7f reduces to a hnear equation in the 
variable x whose solution is 



which in general does not coincide with the distance d between the two pivots. 
Therefore for an ideal Kater pendulum the answers to the previous questions are: 

(1) there are at most three possible positions of the movable weight which determine a 
good mass configuration; 

(2) no; there always exists a good mass configuration whose associated equivalent length 
does not coincide with the distance between pivots. 

An immediate consequence is that a parabolic fitting of the empirical data {x, T) can't 
be the best fitting since two parabolas never meet at three points! Moreover in some 
particular cases a parabolic fitting may cause strong distortions in determining good 
mass configurations. For example: 

• if either d/\/2 < a or (8) is not satisfied, the ideal Kater pendulum admits a unique 
good mass configuration; typically a parabolic fitting of data in this situation gives 
parabolas meeting only at imaginary points and the procedure stops; 

• if a ~ ±d/2 then the first solution of (9) is quite near to one of the two further 
symmetric solutions; a parabolic fitting of data gives only two intersection points 
but we do not know if one (and which one?) of them is nearer to the position 
associated with / than to the one associated with /'; in this situation also / ~ /' but 
they are not equal; then associating I with a so determined good mass configuration 
may cause a relevant error in the final value of g. 

One may object that we are discussing an empiric procedure by means of an ideal pen- 
dulum. In particular the position x — a ior the movable mass gives the completely 
symmetric mass configuration with respect to the middle point of the ideal pendulum 
bar. When a physical pendulum with ruf ^ is considered, what is such a mass 
configuration? Does it occur again? 

The answer is "yes" . The key observation is that, for both pivots, the variable position 
x of the movable mass and the resulting period T of small oscillations are related by a 
cubic expressions of the following type (period-distance relations) 




and the associated equivalent length is 




4a 



ax^ + 6x + c = + dxT'^ , 



(13) 



where the coefficients a, 6, c, d depend on the pendulum parameters. This is precisely 
what Shedd and Birchby pointed out in their papers [2, 3, 4] giving theoretical and 
empirical evidences: they called the two (one for each pivot) equations (13) the equations 
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of the reversible pendulum (see equations (10) and (11) of their first paper). Here we 
will refer to (13) as the cubic period-distance relation of the physical Kater pendulum 
considered. Note that only coefficients a, b, c, d depend on the pendulum parameters, 
while the polynomial type of equation (13) does not depend on the choice of the pendulum. 
Thus, we can reduce the search for good mass configurations to a simple cubic equation 
similar to Eq. (9). 

A first point in the present paper is to give a mathematically rigorous proof of the 
following 

Theorem 1.1. Let pi {x, y) ,p2 {x, y) be the following cubic polynomials 

Pi{x,y) ^ Aix'^ + BiX + Ci-y'^ - Dixy'^ , i = l,2. 

where Ai, Bi,Ci, Di are real coefficients and Di ^ D2 ■ Then they admit always two real 
common roots and two pairs of complex conjugated common roots which may be real under 
suitable conditions on coefficients A^, Bi,Ci, D^. Thinking them as points in the complex 
plane {x, y) they are symmetric three by three with respect to the x-axis. Moreover these 
are all the common roots they can admit (that is: all the further common roots are "at 
infinity"). 

This algebraic result leads to the following physical statement: 

Corollciry 1.2. A physical Kater pendulum, with a "sufficiently long" bar, admits al- 
ways a "good" mass configuration whose associated equivalent length does not in general 
coincide with the distance between the pivots. 

Under "suitable conditions" on the pendulum parameters, it may admit two further 
"good" mass configurations. They correspond to symmetric positions of the movable mass, 
with respect to the middle point of the bar (if also the pivots are symmetrically located). 
They admit a common associated equivalent length which is precisely the distance between 
the pivots. 

Moreover the pendulum can't admit any further good mass configuration. 

We will specify the meaning to the vague expressions "sufficiently long" and "suitable 
conditions" . 

Although Shedd and Birchby knew in practice the content of the previous statement 
(they actually wrote down all the three good mass configurations in period-distance terms 
- see formulas (27) of their first paper) they couldn't give a rigorous proof of it. They 
studied the geometry of the curves determined by the cubic period-distance relations 
by means of an old and non-standard "Newton's classification". They then arrived to 
conclude that (see the bottom lines of p. 281 in their first paper): 
"Of the nine possible intersections of two cubic curves, in the present case three are 
imaginary or at infinity, three belong to the condition that T is negative, and three 
belong to positive values of T, and can hence be experimentally realized. " 
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This conclusion does not exclude that, under some suitable conditions on the pendulum 
parameters, at least two of the three "imaginary or at infinity" intersections may become 
real and maybe physically realizable giving more than three good mass configurations. 
Actually we will see that these three intersections are not imaginary but definitely "at 
infinity" and they can never give physical results. 

A second aim of the present paper is to observe that the best fitting of empirical data 
{x,T) is then given by two cubic curves of type (13) instead of two parabolas of type 
(5). We will support this remark by experimental evidence for a real compound Kater 
pendulum. 

The paper is organised as follows. Section 2 is devoted to prove Theorem 1.1. and 
the physical statement of Corollary 1.2.. Here we set the main notation and describe the 
physics of a real Kater pendulum. The proof of Theorem 1.1. is based on elementary 

elements of complex algebraic and projective geometry. Anyway a non-interested reader 
may skip it without losing any useful element to understand what follows. In Sec. 3 we 
describe an effective experiment. Section 4 is devoted to the analysis of experimental data 
by a linear fit of the period-distance cubics. An estimate of their intersection points gives 
the good mass configurations and the associated periods for our real Kater pendulum. 
Then the value of g. A comparison with a parabolic fitting of data is then given. 

Appendix A is devoted to discuss the "suitable conditions" on the pendulum param- 
eters under which the pendulum admits all the possible good mass configurations (see 
Corollary 1.2.). In Appendix B we collect some further numerical methods to analyse 
our empirical data. 



2 Physics of the Kater reversible pendulum 

Notation. Consider a physical Kater pendulum composed of a rigid bar equipped with 
two weights (see Fig. 2 and Fig. 3). The pendulum can be suspended by two knife-edges, 
ci and C2, symmetrically located on the bar. The weight m/ is placed in a fixed position 
which is not between the knives. The other one, rrim, can be moved along the bar. Small 
oscillations of the pendulum are parameterised by an angle (f such that (f ~ sin (p, that 
is, </7^ ~ 0. The equation of motion of the pendulum is then given by 

^+^^ = 0, (14) 

where g is the earth's apparent gravitational acceleration, m is the total mass of the 
pendulum, hi is the distance of the center of mass from the knife-edge q, and is the 
moment of inertia about q. 

The Steiner's theorem [6] asserts that 

Ii^Io + mh] , (15) 



where /q is the moment of inertia with respect to the center of mass. The associated 
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movable mass 



fixed moss 




fixed knife edge 



fixed knife edge 



Fig. 2 Detailed side view of the Kater pendulum (not to scale). 




Fig. 3 Front view of the Kater pendulum (not to scale). The pendulum swings in the plane of 
the picture; its pivot can be inverted. 

period of small oscillations is 



^ 27r h 



mghi 



(16) 
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Equation (16) implies that the Kater pendulum oscillates with the same period as a 
simple pendulum whose length is given by 

* mhi mhi 

Assume now that the movable mass nim is placed at a point xq on the bar such that 

Ti = ^2 = T{xo) . (18) 

Such a point will be called a characteristic position of the pendulum. Equation (18) can 
be satisfied if and only if I1 — I2 — 1{xq). The length I — 1{xq) will be called the character- 
istic length of the pendulum associated with the characteristic position xq. Analogously 
the associated periods T{xq.) will be the characteristic periods of the pendulum. The 
knowledge of 1{xq.) and T{xq.) for each j = 1, 2, 3 yields the value of g from the relation 

T^2n^, (19) 

and therefore 

47r2 

9 = Y^- (20) 

The variable position of rrim is described by a linear coordinate x having origin at Ci. 
Then C2 is the point x = d > (see Fig. 2) while the fixed weight is placed at Xf such 

that {d- L)/2 <Xf <0. 

The movable and fixed weights are composed of disks whose radii are given respectively 
by Tra and r/. 

L is the length of the pendulum bar. 

Therefore, the distance h between the pendulum center of mass and the origin ci 
depends on the position x of and is given by: 

rub + mf + mm 

where rub is the mass of the bar. Set 

m^mb + mf + mm, (22) 

and 

^^jmb + Xfruf 
m 

Then h can be rewrite as follows 

h^K+'^x. (24) 
m 

The moment of inertia /q is then given by 



/o = (/i - Xff mf + {h- xf mm+ (^h- ^ rUb + Iq , 



(25) 
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Set 



^0 = + + (26) 



(27) 



where 

Ii = IS + m,(x,-Kf + mJi-K] +m„K' 
and Iq can be rewritten as follows: 

h - rn^ "^ ~ ""^ x^ - 2mmKx + . (28) 
m 

From Eq. (17) the condition (18) is satisfied if and only if 

Iq + rnhj ^ Iq + m/i| 
m/ii m/i2 ' 

which is equivalent to requiring that 

{hi - h2) {mhih2 - /o) = . (30) 

Prom Eq. (24) we have 

h.^h^K+'^x (31) 
m 

h^^d-h^d-K-—x, (32) 

m 

and we get the first characteristic position by imposing that hi — h2, that is, 

xo,^U^{d- 2xf) . (33) 
2 2mm 

Two additional characteristic positions can be obtained by the second factor in Eq. (30). 
By letting m/ii/i2 — /q = and expressing Iq as in Eq. (28), we have 

x--dx- "^^'-"^'^ + '" ^0, (34) 



whose solutions are 



d 1 /„ mK^-mdK + Il) 
-o. = 2 + 2r + 4 (35) 



ci 1 , ^ mK^-mdK + I'^ 

''^=2- W ^ ^n. ■ ^''^ 

To determine the associated characteristic lengths I^xq.) use Eqs. (17), (29), and (31). 
It follows that 

7o + m + ^xo^y 
^(^°''' mK + mmXQ. ' ^^^^ 
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It is then easy to observe that l{xo^) and l{xo^) are equal and constant because Xq^ and 
Xog are symmetric. Precisely 

^(^02) = ^(^03) = hi + h2 = d, (38) 

and they do not depend on the other physical parameters of the pendulum. On the 
contrary, this is not true for l{xo^) because 

U^^^ ) = - + 2^ + ^/ i^ni + rrif) {d ~ 2xff ^ 
^ 2 md 2mm'md 

The reader may compare the characteristic positions (33), (35) and the associated char- 
acteristic lengths (39), (38), now obtained, with those given in Eq. (27) by Shedd and 
Birchby[2]. 

Moreover the period-distance relations of the pendulum (what Shedd and Birchby 
called "the equations of the Kater pendulum" [2]) can be obtained by Eq. (29) when hi 
and /i2 are expressed as in Eqs. (31). When the pendulum oscillates about the pivot q, 
the period Ti and the distance x results to be related by the following cubic relations 



where 



Ax^ + BiX + Q = T^ + DixTf , i = 1, 2 , (40) 



^ = ^ (41) 
gmK 

Si = 



and 



A, = (44) 

gm{d-K) ^ ' 

Sn^nimd 

All the possible characteristic positions are then given by all the common roots of Eqs. 
(40). 

Proof of Theorem 1.1.. For more details on the mathematics here involved see, for 
instance, Harris[7] or Shafarevich [8] among other introductory textbooks on algebraic 
geometry. 
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Consider {x,y) as coordinates of points in the complex affine plane C^. Then equa- 
tions pi {x, y) — and p2 {x, y) — give two cubic complex algebraic curves, Ci and C2, 
whose intersection points are precisely the common roots of pi and p2- We can compact- 
ify by "adding a line at infinity": this procedure produces the complex projective 
plane P^. More precisely, we can consider our complex variables x and y to be a ratio 
of further variables, that is, 

X Y 
X = — and y = ^ ■ (48) 

The equations defining Ci and C2 multiphed by become the following 

AiX^Z + BiXZ^ + CiZ^ = Y'^Z + DiXY^ (49) 
A2X^Z + B2XZ^ + C2Z^ = Y'^Z + D2XY^ . (50) 

which are the defining equations of the projective completions Ci and C2, respectively. 
The main ingredient of the present proof is the following 

Theorem 2.1. (Bezout) Given two distinct irreducible complex algebraic plane curves 
of degree di and d2, their projective completions admits a finite number of intersection 
points. Precisely if every intersection point is counted with its algebraic multiplicity then 
this number is did2. 

In particular the projective completions Ci and C2 meet in 9 points, counted with their 
algebraic multiplicities. The Bezout theorem is a consequence of the Fundamental The- 
orem of Algebra which asserts that on the field C of complex numbers every polynomial 
admits as many roots as its degree. 

The first step is to study the intersections "at infinity" , that is, which belong to the added 
"line at infinity." The equation of this line is Z = and by Eq. (49) it intersects both 
our cubics at yoo (that is, the point X = Z = which is the infinity point of the affine 
y-axis x = 0) and in x^ (that is, the point Y = Z = which is the infinity point of the 
affine x-axis y = 0). Both of these are inflection points for Ci and €2- In y^ the inflection 
tangent line of Ci is given by 

ti : DiX + Z = , 

while the inflection tangent line of C2 is 

t2 : D2X + Z = . 

They cannot coincide since Di 7^ D2. Therefore yoo is a simple intersection point of 
our cubics, that is, it admits intersection multiplicity 1. On the other hand, in Xoo both 
Ci and C2 have the same inflection tangent line which is the infinity line Z — 0. Then 
has intersection multiplicity 2. Consequently these infinity points count 3 of the 9 
intersection points. The remaining 6 intersections must be affine, that is, they cannot 
belong to the compactifying line at infinity. 
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To find them note that for i — 1,2 

UnCi^yoo, (51) 

with intersection multiphcity 3 because it is an infiection point for d with tangent fine 
ti- On the other hand 

tinCs = {1/00,^1,^2} (52) 
t2nCi = {y^,Qi,Q2} (53) 

where Ph /=yoo, Qk Az/oo, and Ph /^Qk, because ti and t2 are always distinct. Therefore, 
the affine intersection points of Ci and C2 cannot belong to the lines ti and t2, and they 
can be recovered by studying the common solutions to the following equations 

2 Aix^ + BiX + Ci 

2 _ A2X^ + B2X + C2 

because those points do not make the denominators vanish. So they are reduced to the 
roots of the following cubic equation 

{Aix^ + Bix + Ci)(l + D2X) = {A2X^ + B2X + C2)(l + Dix) . (56) 

It is a cubic equation with real coefficients. Therefore it admits 3 complex roots one 
of which is surely a real number. The remaining two roots are necessarily complex 
conjugated: their reality depends on the coefficients Ai, Bi,Ci, Di. 

Proof of Corollary 1.2.. Recall the cubic period-distance relations (40). Setting 
Ti — T2 = y they are represented by the two cubic curves Ci and Cg whose coefficients 
are assigned by formulas (41) and (44), respectively. Note that they are real numbers 
and Di 7^ D2 since d ^ 0. The hypothesis of Theorem 1.1. are then satisfied and the 
characteristic positions of the pendulum must be represented by the real affine intersec- 
tion points admitting y > 0. To conclude the proof observe that Eq. (56) divided out by 
Aii'^/gmK {d — K) gives exactly the cubic equation (30). The real root is then given by 
(33) and it always occurs when 

'^\d-2xf \ <L 
mm 

The remaining two roots are then assigned by (35). A discussion of their reality is given 
in Appendix A. 

3 The experiment 



The physical parameters characterising our pendulum are given in Table 1; the digits in 
parentheses indicate the uncertainties in the last digit. 
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Table 1 The physical parameters characterising the pcnduhim. 



mm (g) ruf (g) rub (g) 


Xf (cm) 


/ (cm) 


d (cm) 


rf (cm) (cm) 


1399(1) 1006(1) 1249(1) 


-26.73(1) 


167.0(1) 


99.3(1) 


5.11(1) 5.12(1) 



The bar length is measured by means of a ruler whose accuracy is ±1 mm. The radii 
rm and r/ and the position Xf are measured by a Vernier caliper accurate to ±0.1 mm. 
The masses nib, nT-m, and ruf are determined by means of a precision balance accurate 
to one gram. With reference to the structural conditions in Appendix A, we are in the 
case 3.b i.e. all the three characteristic positions occurs and precisely Xo^,Xo^ are placed 
between the knives while xq^ is on the opposite side of the bar with respect to m/. By 
recaUing Eqs. (33) and (35), we expect that 

xoi = (104.57 ±0.11) cm (57) 
= (61.74 ± 0.40) cm (58) 
X03 = (37.56 ± 0.31) cm , (59) 

with associated characteristic lengths 

/ (,xoJ = (121.44 ± 0.09) cm (60) 
/ (xoj = / (X03) = d= (99.3 ± 0.1) cm . (61) 

Throughout the experiment the movable mass m„j will be placed in successive posi- 
tions, generally 10 cm from each other, except near the theoretical characteristic positions 
(57) where the distances decrease (see the second column in Table 2) * . The period of 
small oscillation about the two pivots are measured for all those positions of m^. These 
periods are measured by recording the time of each of 9 consecutive oscillations when 
the pendulum starts from the angle cpo ^ 6° ± 1°. For this purpose we used a photogate 
timed by an electronic digital counter ^ . We repeated the procedure for 18 positions of 
irim, at first with respect to Ci and then C2. The average of the 9 values is taken to be the 
period at the given position of rrim whose error is given by half of its maximum excursion, 
that is, ~ 0.0018 s. The initial angle (po is sufficiently small that an equation similar to 
Eq. (14) is valid. By expanding an elliptic integral in a power series, it is possible to 
approximately express the period associated with the exact equation of pendulum motion 

(fi + ^^y^ sin = (62) 

by adding corrective terms, [9, 6] to the period expression given in Eq. (19). In the next 

section we will evaluate such a correction. The results are reported in Table 2. 

* We did not choose positions too close to the estimated characteristic positions to prevent the casual 

occurrence of coincident period measures about the two pivots. In fact our distance measures are effected 

by an error of « ±1 mm. Such an error would cause a strong distortion in determining the empirical 

characteristic positions. One of them would be directly determined by direct measure and its error would 

not be lessened by the fitting procedure. 

t The resolution of the LEYBOLD-LH model is 0.1 ms. 
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Table 2 The experimental data. 



X (cmj 


Ti (s) 




10 


2.3613 


2.0615 


20 


2.1492 


2.0337 


30 


2.0363 


2.0089 


35 


2.0016 


1.9999 


40 


1.9838 


1.9931 


45 


1.9733 


1.9911 


50 


1.9754 


1.9894 


55 


1.9799 


1.9908 


58 


1.9846 


1.9924 


65 


2.0055 


2.0002 


68 


2.0173 


2.0064 


75 


2.0470 


2.0273 


UK 
oO 


9 nQ'?Q 

z.uyoy 


z.yjo 1 o 


90 


2.1224 


2.0969 


92 


2.1334 


2.1071 


106 


2.2178 


2.2174 


110 


2.2441 


2.2589 


120 


2.3078 


2.3776 



4 The linear fitting procedure 

We now describe a linear fitting procedure used to fit the experimental data listed in 
Table 2 and empirically determine the characteristic positions. The numerical compu- 
tations were obtained using MAPLE ^ and some FORTRAN code ^ . From a numerical 
point of view we should fit the data by cubic polynomials like those in Eq. (40). Such a 
fitting can be treated linearly because the coefficients Di and D2 may be determined a 
priori by Eqs. (41) and (44) which involve only the known physical parameters listed in 
Table 1. We obtain 

Di = (3.983 ± 0.01) 10-^ cm"^ (63) 
D2 = (-4.2689 ± 0.0047) 10"^ cm"^ . (64) 
^ We used Maple V, Release 5.1 by Waterloo Maple Inc. 

^ The FORTRAN codes employed subroutines from Rcf. [14] and the numerical package NAG-Mark 14. 
The plotting package is PGPLOT 5.2 developed by T. J. Pearson. 
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Table 3 Coefficients of the cubic curve Ci estimated by the linear method. 
Ai (0.001607 ± 0.000003) cm,-^ 
Bi cm~^ 

Ci (7.641 ± 0.011) s2 

Table 4 Coefficients of the cubic curve C2 estimated by the linear method. 
A2 (0.000172 ± 0.000002) cm-^ 
B2 (-0.03422 ± 0.00031 ) cm'^ 
C2 (4.393 ±0.01) 



Table 5 x'^ and critical p vahw^s for hncar fitting by cn})ic curves 



Ci, degrees of freedom = 16 


; C2 degrees of freedom = 15 


2.69 


1.57 


J^\'^{x,15)dx 0.00008 


0.000005 



We can obtain the desired fitting of the data obtained in Sec. 3 by applying the least 
squares method to the following function 

18 / rp2 _ Aix\+BiXh+C^ \ ^ 

where (xh,Th^i) are the data of the ith set in Table 2 H . 

Two sources of error with period measurements were considered 

• In any position and for both pivots, we considered the standard deviation of the 9 
period electronic measurements, varying from 0.0003 s to 0.0036 s . 

• The systematic error that formula 70 introduces on the data. For example when 
T=2.3s (the maximum period here analysed) and (po ~ 12° the shift introduced on 
T is 0.006 s . 

After this analysis we considered = 0.006 s as estimated error for period measurements. 

The obtained results are reported in Table 3 , Table 4 and visualised in Fig. 4. 

The merit function the associated p-values are reported in Table 5 and each 

of them have to be understood as the maximum probability to obtain a better fitting. 

The estimated cubic coefficients of Table 3 and Table 4 allow us to evaluate the 
characteristic positions and the associated characteristic periods by intersecting their 

II We also know that Bi =0 and the number of coefficients to be estimated by the fitting procedure can 
be reduced. 



M. Rosrii and L. ZaiiiiKUti /' Cciilral Eiiroijcaii S(:i(>iic(> Journals 1 (2003) 1 23 



17 



(0 




X [cm] 



Fig. 4 Theoretical cubics (dotted line), fitted cubics (full line), and experimental data (filled 
points) . The experimental errors are much smaller than the filled points drawn, so they are not 
visible within this plot. 



Table 6 Esl iniau d intersection points of fitting cubic curves. 



(xoi , r(xoi )) (106.015 cm, 2.2184 s) 
(xo2,r(xo2)) (62.541 cm, 1.9973 s) 
(xo3,T(xo3)) (35.779 cm, 1.9998 s) 



upper branches ** We obtain a cubic equation whose numerical solutions are reported in 
Table 6. 

Refer to Eqs. (20) and (60) to compute the associated values of g. We have 



9i 



92 = 4:71 



T(xoJ2 

2 ^C^O^) 



(66) 
(67) 



** These two cubic curves represent the period-distance relations in the plane {x, T) when oscillations 
are considered about ci or C2 respectively. Then their common points coordinates give the characteristic 
positions of the pendulum and the associated periods. We have already observed in Sec. 2 that these 
cubics are symmetrical with respect to the x-axis. More precisely each of them is composed by two 
symmetrical branches. The branches lying under the x-axis are not physically interesting since their 
period coordinate T is negative. Therefore the only interesting common points of these two cubic curves 
are the intersection points of their upper branches. 
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Table 7 Values of g obtained by formulas (66) and (69). 



gi (974.15 ±2.72) cms^^ 

g2 (982.65 it 3.11) cms-^ 

gs (980.20 ±3.1) cms-^ 

g (979.00 ± 1.72) cms-^ 



and their numerical values are listed in Table 7. 
Their average gives 

g = (979.00 ± 1.72) cm s"^ . (69) 

where the uncertainty is found implementing the error propagation equation (often called 
law of errors of Gauss) when the covariant terms are neglected (see equation (3.14) in [11]). 
We now consider the correction arising from the approximation of the exact equation 
of pendulum motion (62) already mentioned at the end of the previous section. This 
correction gives: [9, 6] 



T = 27ry^(l + l/16(^o') , 



(70) 



and 

g = 47:'^ (l + 1/16 ipo'f . (71) 

A small increase in the value of g is evident from Eq. (71), and we will refer to it as the 
finite amplitude correction (f.a.c). 

With the data listed in (7) we obtain 

g^ = (980.34 ± 1.74) cms~^ , (72) 

which is the gravity acceleration increased by the f.a.c. An accurate measure of the value 
of g in Turin [10] gives 

g^ = 980.534099(4) cm s'^ . (73) 

This value will be considered as the "true value" of the acceleration due to the earth's 
apparent gravity field in Turin . By comparing it with 'g_^, we see that our measurement 
is — 191ppm smaller than the "true value." 

Note that the considered Kater pendulum admits characteristic positions sufficiently 
distant from each other (see formulas (57) and data collected in Table 6). Then it can be 

Further references for accurate measurements of g are the following. A world-wide survey of all the 
apparent gravity measurements (see <http://bgi.cnes.fr>) gives for Turin g = 980.5495 cm s~^; this 
value differs from by 16ppm. An analytical formula provided by the U.S. Geological Survey [12] 
needs two input parameters, the local height above sea level and latitude, which in our case are 236 m 
and 45.05333° respectively to give the local value of g. For Turin this formula gives g = 980.5937cm s~^, 
which differs from 51T by 61 ppm. 
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Table 8 Coefficients of the first fitting parabola Vi . 

Ai (0.000180 ± 0.000002) s cwt'^ 
Bi (-0.01959 ± 0.00017) s cm,-^ 
Ci (2.494 ± 0.003) s 



Table 9 Coefficients of the second fitting parabola V2- 

A2 (0.000054 ± 0.000002) s cmr'^ 
B2 (-0.00517 ± 0.00017) s cmT^ 
C2 (2.113 ± 0.004) .s 



Table 10 and critical p- values for linear fitting by parabolas , 
first 13 data , degrees of freedom =10 

Vi V2 

893 12.5 
Sfx\x,lh)dx 1 0.74 



considered sufficiently "well-assembled", which means that a parabolic fitting (of type 
(5)), of the empirical data {xh, T^^i) collected in Table 2, should give a sufficiently precise 
evaluation of characteristic positions xq^^xq^. As before we apply the least square method 
to the following function , but the sum is now extended to the first 13 entries of Table 2 
in order to exclude the first intersection 

h=l \ (^T J 

The coefficients of the fitting parabolas are reported in Tables( 8) and ( 9) ; the ^^iid 
the associated p-values are reported in Table 10 , in comparison the parabohc fit gives 
very bad results. . 

Their intersection points are given in Table (11). 

We get then the following two evaluations of g: 

g' = 47r^ ^, ^ ,^ = (960.08 ± 3.00) cm s"^ (75) 



T(xo. 



'2 



20 



M. Rossi and L. Zaiiiiictli / (\>iitral Euroi)caii Science Journals 1 (20U3) 1 25 



Table 11 Estimated intersection points of fitting parabolas. 

{xo^,T{xo^)) (72.296cm : 2.0207s) 

(,To;,,.T(;Tn,)) (41.709r7?? : 1.9908.s) 



/ = 47r^ ' = (989.11 ± 3.14) cms-^ . (76) 

Their average gives 

9paraboiic = (974.60 ± 2.18) . (77) 

A comparison with g in (69) and gx in (73), gives a clear evidence of the better efficiency 
of a cubic fit with respect to a parabohc one. 



5 Conclusions 

We summarise the main results of our theoretical and numerical analysis 

(1) The three solutions of the Katcr pendulum concerning the distance-period relation- 
ship discovered in 1907 by Shedd and Birchby in 1907 [2, 3, 4] are classified in a 
modern context 

(2) The first solution of the distance-period relationship allows to deduce a new for- 
mula of g via the second equivalent length both in the idealized pendulum and in a 
commercial Kater pendulum, see respectively formula(12) and (39) 

(3) One of the main targets of our work "the evaluation of g" gives oscillating results 

• : our best numerical fit to produces (the linear fit -|- non-linear correction) 
a value of g that is 191 ppm smaller than the "true vale" 

• : our worst fit to (the non-linear fit + non-linear correction) gives a value 
of g that is 1978 ppm smaller than the "true vale" 

(4) Concerning the fit to T through a parabola we obtain high values of (x^=893 for 
Ci and x^=12.51 for C2) with respect to the hnear fit to (x^=2.69 for Ci and 
X^=1.57 for C2). These high values of allow to rule out the physical significance 
of this type of fit. 
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A Reality of characteristic positions and structural conditions 



The characteristic positions of our pendulum are given by Eqs. (33) and (35). The former, 
, is always real. On the other hand Xq^ and Xq^ are real if and only if the square roots 
in Eq. (35) are real i.e.e if and only if 



rumd^ + AmdK - AmK^ - 4/n > 



x'j — dxf — 



(mm + mi,)d^ - AI'^ 
Aruf 



< 0. 



(A.1) 



The latter are the "suitable conditions" on the pendulum parameters of Corollary 1.2.. 

To avoid the overlapping of rrif with ci we have to impose that Xf < —rf. Then 
Eq. (A.l) is equivalent to requiring that 



md'^ - AI'' > and 



1 md2-4K 



mf 



< Xf < — rf 



(A.2) 



Note that the condition on the right in Eq. (A.2) is not empty if 

d 1 md2 - 47/,' 



2 2\ ruf 



< - 



d>2 



mfTf 



In particular, the latter ensures that the left condition in Eq. (A.2) is also satisfied because 

-'o 



2 ^ mfr'j + Jq 



(A.3) 



md^ - 4/^' > 



d>2\ 



m 



(A.4) 



To avoid the overlapping of rrim with the knife-edges it follows that either < a; < d—rr, 
OT d + rm ^ X < After some algebra we get the following results. 



Assume that > n^-f and set 

.rrimrm + mfTf 



Ml 

Si 
S2 



mm - mf 

ylTlb + mm 

ruf - mm 



I 



+ 



\\m\) + mm) TJib + m^ 



2m/ 

d_l 
2~2\ 

d 1 



, , mm 

d H Tr, 

mf 



md^ - Aid 
mf 



mfd^ + mbd^ + Ammrmd - AmmTm - 4/o 



^ 2 2\ 
Then we have the following possibilities: 



mf 



(A.5) 
(A.6) 
(A.7) 
(A.8) 

(A.9) 
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(1) d < min(Mi, M2): in this case the pendulum admits only one characteristic position 
given by xq^ because by Eqs. (A. 2) and (A. 3), xq^ and X03 are not real; xq^ is not between 
the knives, but occurs on the opposite side of the bar with respect to m/; the system is 
in an almost symmetrical mass configuration of the pendulum. 

(2) For min(Mi, M2) < d < max(Mi, M2), we have the following possibilities: 

(2a) If Ml < M2: we get only the characteristic position a^Oi which is between the 
knives if and only if Si < Xf < —rf \ otherwise, x^^ is placed like in (1). 

(2b) If M2 < Ml, we get all the characteristic positions; is like in (1) and a;o2,a;o3 
occur between the knives if and only ii S2 < Xf < min(S'3, — r/). 

(3) For max(Mi , M2) < d, the pendulum admits all the characteristic positions xq^ , xq^ , X03 
which are placed as follows: 

(3a) Only xq^ is placed between the knives when either 

Si < S2 and Si<Xf<S2, (A.IO) 

or 

5*3 < —rf and max (5*1, ^3) < Xf < —Vf. (A. 11) 

(In particular, if Ss < Si, we can also assume the position Xf — Si for the fixed weight 
ruf.) 

(3b) only placed between the knives when S2 < Si and S2 < xj < Si. 

(3c) We obtain all the possible characteristic positions Xq^ , Xq^ , Xq^ between the knives 
when 

max(S'i, 5*2) < Xf < min(S'3, —rf) . (A. 12) 

In the concrete case considered in the Section 3 we have 



Ml 


= (62.61 ± 0.25) cm 


(A.13) 


M2 


= (70.87 ± 0.05) cm 


(A.14) 


Si 


= (-12.28 ±0.08) cm 


(A.15) 


S2 


= (-28.049 ± 0.075) cm 


(A.16) 


Ss 


= (-7.62 ±0.06) cm. 


(A.17) 



where the uncertainty is found by applying the law of errors of Gauss with the uncertain- 
ties listed in Table 1. Therefore we are in the case 3.b. 

B Further numerical methods 

We outline three additional numerical methods that may be applied to analyse experimen- 
tal data. The final results obtained by means of each method are reported in Table B.l. 

B.l The non-linear method 

In the fitting procedure of data reported in Table 2, all the coefficients Ai,Bi, Ci and Di 
are considered as unknown parameters to be estimated. Therefore a fitting procedure 
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Table B.l Average values g and corrected values g_^_ (by f.a.c). 



algorithm 


g 


54- 


linear fitting by parabolas 
linear fitting by cubics 
non-linear fit 
Cramer interpolation 
Spline interpolation 


{974.6 ±2.17) cm s-^ 
(979.00 ± 1.72) cm s-2 
(977.25 ± 1.71) cm s-2 
(980.06 ± 4.88) cms'^ 
(979.52 ± 1.73) cm 


(975.93 ±2.20) cm s-2 
(980.34 ± 1.74) cm 
(978.25 ± 1.74) cm 
(981.40 ± 4.89) cms-2 
(980.86 ± 1.74) cms-2 



performed by means of cubic polynomials like those in Eq. (40) is necessarily a non-linear 
one. We want to apply the least square method to minimise the following functions 

X, (A, A, a, A) = E - , (B.l) 

which are non-linear in the unknown coefficients. The procedure is to apply the NAG- 
Markl4 subroutine E04FDF to find an unconstrained minimum of a sum of 18 nonlinear 
functions in 4 variables (see Ref. [13]). 

The final value of g is reported in the Table B.l. 

B.2 The Cramer interpolation method 

We present here a method that reduces our analysis in a local neighbourhood of the 
estimated characteristic positions where a cubic behaviour of the fitting curves is imposed. 

Prom Eqs. (41), (44), and (63) we know that D\ and D2 are completely determined 
by the pendulum parameters. Moreover, from Eqs. (41) we know that — 0. So to 
recover the remaining coefficients of C\ and C2, we need to interpolate two points of the 
ffist set of data in Table 2 and three points of the second one respectively. We have to 
solve a 2 X 2 and a 3 x 3 linear system by applying the Cramer theorem (which is the 
most practical method for solving a linear system of equations) . If we choose data points 
that are close to a characteristic position, then the nearest point in C\ fl C2 to the chosen 
data will give an empirical estimation of such a characteristic position and its associated 
period. An iterated application of this procedure will produce a distribution of periods 
and we may obtain g from the mean value and its statistical error from the standard 
deviation (see the last line of Table B.2). 

The results obtained for every interpolation are reported in Table B.2; the chosen 
data points in the second and third columns are enumerated as they appear in Table 2. 

B.3 Cubic Sphne Interpolation 

The last data analysis method to be proposed is the cubic spline interpolation (subrou- 
tine SPLINE and SPLINT from Numerical Recipes II). Once the three intersections are 



24 



M. Rossi and L. Zaiiiiictli / (\>iitral Euroi)caii Science Journals 1 (2003) 1 25 



Table B.2 The Cramer interpolation method. 



Cramer Method 



Chctr. 


Chosen data 


Intersection 






Char. 


0" 

6 








oei ICb 1 


oeiico z 


JT UbiliUli 


Period 




1 PTi P'f, n 








1fe\ 


10,1(5 


10,1 / ,1(5 


1 77^ nm 
lUO. MO CIIl 


+11 — 

Ll,l — 


9 91 7 


c 


1 91 44 rm 


Q7^ 7S 


i^iii 


-2 




10, 1 ( 


10,1D,1 1 


lUD.oDU cni 


+ 1 

L1,Z — 


9 991 


c 


1 91 44 PTn 


Q71 Q9 




-2 


irr\ 


ID, io 


1 K-1 7-1 8 
ID, 1 ( , lo 


lUD.iuo cm 


ti,o — 


9 918 
z.zio 


c 


TZiT.4:4: t^ill 


Q74 ns 


r*vn Q 


-2 


XQi 


10, lo 


10, ID, lo 


lUD.loy cm 


ft A— 


9 91 Q 
z.ziy 


Q 


1 91 44 cm 


Q73 44 


r'TTi c:~ 
i^iii 


-2 




7- Q 


7- 8- Q 


fi9 7SQ /^m 

oz. ( oy cm 


t2,l= 


1.996 


C 


QQ '^n rm 


Q83 77 




-2 


1fe\ 


8;10 


8; 9;10 


62.056 cm 


t2,2= 


1.996 


C 




Q83 7Q 




-2 




9;11 


9;10;11 


61.962 cm 


t2,3= 


1.996 


s 


99.30 cm 


984.28 


cms" 


-2 


X02 


10;12 


10;11;12 


61.990 cm 


t2,4= 


1.996 


s 


99.30 cm 


984.44 


cms" 


-2 


XOg 


2; 4 


2; 3; 4 


35.477 cm 


t3,l= 


1.999 


s 


99.30 cm 


980.87 


cm s~ 


-2 


XO3 


3; 5 


3; 4; 5 


36.207 cm 


t3,2= 


1.998 


s 


99.30 cm 


981.97 


cm s" 


-2 


XO3 


4; 6 


4; 5; 6 


35.557 cm 


t3,3= 


1.999 


s 


99.30 cm 


981.11 


cm s" 


-2 




5; 7 


5; 6; 7 


36.668 cm 


t3,4= 


1.995 


s 


99.30 cm 


985.38 


cms" 


-2 



5=( 980.06 ± 4.89) cms-^ 



obtained, the procedure was similar to the hnear/nonhnear case and the final value of g 
is reported in Table B.l. 
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